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Abstract: Recent experimental and theoretical determinations of (g — 2)^, QcDhih 2 and 
BF(b — ► 57) place exceedingly tight constraints on the minimal supergravity (mSUGRA) 
model. We advocate relaxing the generational universality of mSUGRA, so that GUT 
scale third generation scalar masses are greater than the (degenerate) first and second 
generation scalar masses (a normal scalar mass hierarchy (NMH)). The non-degeneracy 
allows for a reconciliation of all the above constraints, and also respects FCNC limits from 
Bd — B<i mixing and b — > 57. The NMH SUGRA model leads to the prediction of relatively 
light first and second generation sleptons. This yields large rates for multilepton collider 
signatures at the CERN LHC and also possibly at the Fermilab Tevatron. The spectrum 
of light sleptons should be accessible to a y/s = 0.5 — 1 TeV linear e + e~ collider. 
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1. Introduction 



The minimal supergravity model (mSUGRA) has long served as the paradigm supersym- 
metric model for testing the phenomenological consequences of weak scale supersymmetry ]l[]. 
The model is constructed by starting with the general supergravity Lagrangian of N = 1 
supersymmetric gauge theories^]. The gauge symmetry and matter content of the Stan- 
dard Model are assumed, where now matter fermions comprise just one component of left 
chiral superfields, and the Higgs sector must be expanded by an additional doublet to can- 
cel triangle anomalies. Supersymmetry breaking occurs within a postulated hidden sector, 
and can lead to weak scale soft SUSY breaking terms of mass m = A 2 /Mpi if the SUSY 
breaking parameter A ~ 10 11 GeV. Motivated by the degeneracy solution to the SUSY 
flavor and CP problems^, the simple choice of a fiat Kahler metric can be made. This 
leads at tree level to degenerate scalar masses at some high mass scale, usually taken to 
be Mqut ~ 2 x 10 16 GeV. A simple choice of the gauge kinetic function leads also to 
gaugino mass unification. Taking the limit as Mpi — ► oo while keeping m 3j / 2 fixed leads 
to a globally supersymmetric renormalizable gauge theory where supersymmetry is broken 
by the presence of soft SUSY breaking terms (the MSSM). The weak scale parameters are 
related to those of the GUT scale by RG evolution, and electroweak symmetry is broken 
radiatively. The weak scale sparticle masses and mixings are determined by the parameters 

m 0! mi/2, A , tan/?, sign(fi), (1.1) 

where m,Q is the common scalar mass, is the common gaugino mass, and Aq is the 

common trilinear soft SUSY breaking term all valid at Q = Mqut- Throughout this work, 
we will adopt the value rat = 175 GeV. 

The mSUGRA model as defined above determines the large number of MSSM param- 
eters in terms of just a few parameters stipulated at Q = Mqut, and therefore provides 
a very rich phenomenology in terms of just a few free parameters. However, the model 
has long been susceptible to criticism owing primarily to the assumption of universality, 
especially amongst scalar masses[||]. First, there is no theoretical reason for choosing a 
flat Kahler metric, and in general the scalar masses can be arbitrary. The arbitrary soft 
breaking masses in general lead to violation of constraints from FCNC processes ||, an d 
complex soft breaking parameters can lead to CP violating phenomena at large rates 0. It 
has been shown that even if universality is imposed at tree level, the universality will be 
broken by loop corrections 0]. We note here that there do exist symmetries which guar- 
antee universality of soft breaking terms within a generation, for instance, 50(10) gauge 
symmetry. But then the non-universality would still be manifest as a splitting amongst 
the generations, and with the Higgs boson soft masses. 

In recent years, supersymmetric models have become increasingly constrained by a 
variety of measurements ||. These include determination of the branching fraction BF(b — > 
s l) HE0]> — ' H' El' —I' ^ ne m uon anomalous magnetic moment a M = (g— 2)^/2 fL5|l and the 



relic density of cold dark matter in the universe. The most recent development comes from 
the Muon g — 2 Collaboration, which published results on (g — 2)^ for the negative muon 
along with earlier results on the positive muon. In addition, theoretical determinations of 



- 1 - 



(g — 2) a have been presented by Davier et ai.|l£] and Hagiwara et al. p7|1 which use recent 



data on e + e~ — ► hadrons at low energy to determine the hadronic vacuum polarization 
contribution to the muon magnetic moment. Combining the latest experiment and theory 
numbers, we find the deviation of to be: 



Aa M = (27.1 ± 9.4) x 1(T 10 (Davier et al.) (1.2) 
Aa M = (31.7 ± 9.5) x 10~ 10 (Hagiwara et al). (1.3) 

The Davier ei a/, group also presents a number using r decay data to determine the 
hadronic vacuum polarization, which gives Aa^ = (12.4±8.3) x 10 -10 , i.e. nearly consistent 
with the SM prediction. However, there seems to be growing consensus that the numbers 
using the e + e~ data are more to be trusted, since they offer a direct determination of 
the hadronic vacuum polarization. The ~ 3a deviation in using the e + e~ data can 
be explained in a supersymmetric context if second generation sleptons (smuons and mu 



sneutrinos) and charginos and neutralinos are relatively light [18]. 

In SUSY models with minimal flavor violation, the flavor changing decay b — > 57 can 
still occur. The decay proceeds via Feynman graphs including t\ -jfflx 2 and tH + loops, in 
addition to the SM contribution from a tW loop. The branching fraction BF(b — ► 57) 



has recently been measured by the BELLEQ, CLEO|l^] and ALEPH|TlJ] collaborations. 
Combining statistical and systematic errors in quadrature, these measurements give (3.36± 
0.67) x 10" 4 (BELLE), (3.21 ± 0.51) x 10~ 4 (CLEO) and (3.11 ± 1.07) x 10~ 4 (ALEPH). 
A weighted averaging of these results yields BF(b — > 57) = (3.25 db 0.37) x 10~ 4 . To this 
we should add uncertainty in the theoretical evaluation, which within the SM dominantly 
comes from the scale uncertainty, and is about 12%. 1 Experiment and theory together 
imply the bound, 

BF(b^ 57) = (3.25 ±0.54) x 10" 4 . (1.4) 
The calculation of SM and supersymmetric contributions to BF(b — * sj) used here is based 



upon the program of Ref. [14]. Since the SM value of BF(b — > 57) is relatively close to the 
central measured value, rather large values of 2 , 2 and are expected, which 

help to suppress the SUSY loop contributions to the decay amplitude. 

Finally, the recent precision mapping of anisotropics in the cosmic microwave back- 
ground radiation by the WMAP collaboration have fitted the relic density of cold dark 



matter (CDM) in a ACDM universe[19, 20|. The constraint on the relic density of neu- 



tralinos produced in the early universe is then 

Q^h 2 = 0.1126 ± 0.0081. (1.5) 

The new relic density constraint is highly restrictive. In the mSUGRA model, only a few 
regions give rise to this narrow band of h 2 : 



1 We caution the reader that the SUSY contribution may have a larger theoretical uncertainty, particu- 
larly if tan p is large. An additional theoretical uncertainty that may increase the branching ratio in the 
SM is pointed out in Ref. Ftl. 
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• The bulk region at low mo and low my2, where neutralino annihilation in the early 
universe proceeds dominantly via i-channel slepton exchange. This region is now 
largely excluded (save where it overlaps with the stau co-annihilation region) because 
the WMAP allowed region gives rise to low values of in conflict with LEP2 bounds, 
and also large deviations in the BF(b — > S7). 

• The stau co-annihilation region at low mo where — , so that Z\ — t\ and 

— 7j~ co-annihilation provides a ready sink for neutralinos in the early universe ||21||, 

• The hyperbolic branch/focus point region (HB/FP) at large mo near the boundary of 
parameter space where /j, 2 becomes small[^2|. The Z\ becomes increasingly higgsino- 
like in this region, facilitating its annihilation rate into WW, ZZ and Zh in the early 
universe. 

• The A-resonance annihilation region at large tan/3, where 2m^ ~ m^, m#. In this 
region, the A and H widths can be quite broad (~ 10 — 40 GeV), giving rise to a 
relatively broad A- annihilation "funnel" [^] . 

In addition, other regions can occur such as the h resonance annihilation region (a very 
narrow strip at low rrtx/2 where 2m^ ~ m^) and stop co-annihilation region (at the edge 
of parameter space for very particular values of Aq parameter). 

A recent x 2 analysis of fi^/i 2 , Aa M and BF(b — > sj) has been made[p4|l, and has 
found the HB/FP, A-funnel and stau regions to coincide with all experimental constraints. 
The analysis of Ref. |24j] pre-dated the new experimental and theoretical determinations of 
Aa^ 2 . The improved Aa M numbers, as determined experimentally by E821 and theoreti- 
cally by Davier et al. and by Hagiwara et al, point to a larger deviation from the SM, and 
in particular, to relatively light second generation slepton masses, which do not occur in 
the HB/FP region, or in the large m^j part of the A- funnel or stau co-annihilation region 
of the mSUGRA model. 

Our goal in this paper is first to present updated x 2 fits to the mSUGRA model, 
including the new Aa M determination. While the new Aa M numbers prefer lighter sleptons, 
the BF(b — * 57) numbers prefer heavy squarks, typically in the several TeV range. The 
tension between these two constraints, coupled with the tight limits from WMAP, 

leave only small viable portions of mSUGRA parameter space. We present these results in 
Sec. 

In Sec. IH, we advocate relaxing the universality assumption of the mSUGRA model. 
In this case, constraints from FCNC processes must be addressed. The FCNC constraints 
apply most strongly to first and second generation scalar masses. We explore a scenario in 
which first and second generation scalars remain degenerate, while allowing for a significant 
splitting with third generation scalars. In this case, heavy (multi-TeV) third generation 
scalars are preferred by BF(b — * 57) constraints, while rather light first and second gener- 
ation scalars are preferred by Aa^. The scenario is called the normal scalar mass hierarchy 

2 Explictly, the theory calculation of Narison^HJ using e + e~ — > hadrons data was used, which gave 
Aa,, = (24.1 ± 14) x 1CT 10 , i.e just a 1.7a deviation from the SM prediction. 
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(NMH), to distinguish it from earlier studies which advocated an inverted scalar mass hi- 
erarchy (IMH) p6[| . The scenario respects constraints from Bj — mixing. The resulting 
model parameter space only allows low values of mi^, save for a small portion of the 
HB/FP region. Significant regions of parameter space survive all three constraints as given 
above. Light first and second generation sleptons of mass just a few hundred GeV are 
characteristic to this scenario. 

In Sec. we examine consequences of the NMH scenario for collider experiments. 
When first and second generation sleptons become light, then chargino and neutralino de- 
cays to electrons and muons are enhanced, leading to large rates for distinctive multilepton 
plus jets plus events at the CERN LHC, and possibly observable trilepton rates at 
the Fermilab Tevatron. At a linear e + e~ collider operating at yfs = 0.5 — 1 TeV, first 
and second generation sleptons should be accessible to discovery, while tau sleptons would 
likely be beyond reach. 

We summarize our results in Sec. ||. 



2. x 2 At of mSUGRA model 

Our first objective is to calculate the \ 2 quantity formed from Ao^, BF(b — ► S7) and 
h , as given in Eq'ns pP ). For % h 2 , we only use the upper bound contribution 
to the x 2 ) since only the upper bound from WMAP is strictly applicable, i.e. there may be 
other forms of dark matter in the universe. In addition, we use the SM determination as 
given by Hagiwara et al, which includes the latest fit to the hadronic vacuum polarization 
from e + e~ — ► hadrons at low energy. We generate the mSUGRA model particle spectrum 
using Isajet 7.69|p7|1, which includes two loop RGE running of all couplings and soft SUSY 
breaking terms, minimizes the one loop effective potential at an optimized scale choice 
(which accounts for leading two loop terms), and which implements the complete set of log 
and finite corrections to all sparticle masses 3 . 

Once the superparticle mass spectrum and mixings are generated, we calculate the 
quantity Qgji 2 using the Isared program, which includes relativistic thermal averaging^] 
of all relevant neutralino annihilation and co-annihilation processes as calculated exactly 



at tree level using the CompHEP program[3C]. We evaluate Aa^ as in Ref. and the 



BF(b — > 57) as outlined in Ref. [14], which uses a tower of effective theories approach [plfl 
and includes NLO corrections to the running of Wilson coefficients between Q = M\y and 
Q = m.6|3^], and also includes the full set of NLO QCD corrections to the b — > 57 process 
at the scale Q = m^Q. In addition, in our calculations, we implement the running value 
of the b Yukawa coupling including MSSM threshold corrections for calculations above the 



scale Q = Myy pi]. Our BF(b — > 57) results then agree well with those of Ref. p| at 



both small and large tan f3. 

For every point in mSUGRA model parameter space, we calculate the x 2 quantity 
formed from Aa^, BF(b — ► sj) and fi^/i 2 . For convenience, we plot in the tuq vs. m\/2 
plane the quantity \/x^, which is color coded according to the legend shown in Fig. Ill: i.e. 



3 Good agreement is found between Isajet and programs such as Suspect, SoftSUSY and Spheno in most 
of the mSUGRA model parameter space pq] 
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green regions have relatively low x 2 , while red has a high x 2 an d yellow is intermediate. 
For variations in the two degrees of freedom (mo and mj/j), x 2 = 2.3 (\/x 2 = 1-52) 
corresponds to a la deviation from the combined measurements, x 2 = 6.18 ^\fx 2 = 2.49) 
corresponds to a 2a deviation and x 2 = H-83 (a/x 2 = 3.44) corresponds to a 3a deviation. 

Our first results are shown in Fig. [I], where we plot for Aq = 0, fj, > and tan (3 values 
of 10, 30, 55 and 58. For fi < 0, x 2 fits are almost always worse, since Aa^ favors fi > 0. 
Variations in the Aq parameter typically lead to no quantitative changes, save for those 
special Aq values which drive the top squark mass to very low values. The gray regions are 
excluded theoretically, while the blue regions are excluded by LEP2 searches for chargino 
pair production. In Fig. ||, we present contour plots for BF(b — > 57), a M and S7/i 2 for the 
same parameters as Fig. [l]; these figures complement those of Fig. [l] and explain the details 
of the x 2 behavior. 

In the tan/3 = 10 case, it can be seen that almost all the m,Q vs. m^ plane has very 
large x 2 ■ This arises because in general an overabundance of dark matter is produced in the 
early universe, and the relic density ^^fo 2 is beyond WMAP limits. There is a very narrow 
sliver of yellow at m^ 2 ~ 150 GeV (just beyond the LEP2 limit) where 2m^ — m^, and 
neutralinos can annihilate through the narrow light Higgs resonance. In addition, there is 
an orange/yellow region at high mo at the edge of parameter space (the HB/FP region), 
with an intermediate value of x 2 ■ I n an earlier study [24], this region was found to have a 
low x 2 value. In this case, however, the 3a deviation from the SM of tends to disfavor 
the HB/FP region. In the HB/FP region, sleptons are so heavy (typically 3-5 TeV), that 
SUSY contributions to a M are tiny, and the prediction is that a M should be in near accord 
with the SM calculation. The remaining green region is the narrow sliver that constitutes 
the stau co-annihilation region, barely visible at the left hand edge of parameter space 
adjacent to where f\ becomes the LSP. The overall situation is similar when we move to 
tan (5 = 30: the only low x 2 region is in the stau co-annihilation corridor. 

Once we move to very large tan (3 values, as shown in the third frame, then the A- 
annihilation funnel becomes visible, and some large regions of moderately low x 2 appear 
around mo, m^j ~ 500, 600 GeV and also at 1500, 200 GeV. While the ^-annihilation 
funnel extends over a broad region of parameter space, the upper and lower ends of the 
funnel are disfavored: basically, if sparticles become too heavy (the upper end), then Aa M 
becomes too small, while if sparticles become too light (the lower end), then BF{b — > 57) 
deviates too much from its central value. The final frame shows results for tan j3 = 58, close 
to the point where parameter space begins to collapse due to an inappropriate breakdown 
of electroweak symmetry. The relic density is low over much of the parameter space since 
2m^ ~ rriA, so a region of relatively low x 2 appears at modest mo and myz values. 

Our conclusion for the mSUGRA model is that almost all of parameter space is ex- 
cluded or at least disfavored by the combination of the WMAP Slj^fr 2 limit, the new Aa^ 
value, and the BF(b — > 57) value. The Sl^ft- 2 constraint only allows the several regions 
of parameter space mentioned in Sec. |], while BF(b — > S7) favors large third generation 
squark masses to suppress SUSY contributions to b — > S7 decay, and Aa M favors relatively 
light second generation slepton masses, to give a significant deviation of {g — 2)^ from 
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mSUGRA, tanp=1 0, \i>0, A D =0, m top =1 75 GeV mSUGRA, tanp=30, \i>0, A D =0, m top =1 75 GeV 

e + e input for 8a • LEP2 excluded e + e input for Sa u • LEP2 excluded 
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mSUGRA, tanp=55, fj>0, A Q =0, m top =1 75 GeV mSUGRA, tanp=58, \i>0, A Q =0, m top =1 75 GeV 
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Figure 1: Plot of regions of VX 2 i n the mSUGRA model for Aq — 0, fi > 0, and tan j3 = 10, 30, 55 
and 58. The green regions have low x 2 , while red regions have high \ 2 ■ Yellow is intermediate. 



the SM value. The only surviving regions with relatively low \J~x^ ~ 2 are the stau co- 
annihilation region, and intermediate portions of the ^4-annihilation funnel at very large 
values of tan/3. The previously favored HB/FP region from Ref. [24] is now disfavored 
because it is unable to generate a significant deviation of (g — 2)^ from the SM prediction 
which uses the e + e~ — > hadrons data for evaluating the hadronic vacuum polarization 
graphs. (We note that should the hadronic vacuum polarization determination using r 
decay data turn out to be correct, then the HB/FP region will appear in a more favorable 
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Figure 2: Contour levels for BF(b — > sj), a M and Slh 2 in the mSUGRA model for A® = 0, [i > 0, 
and tan/3 = 10, 30, 55 and 58. Magenta region excluded by LEP2 searches. 

light!) 

3. x 2 At of NMH SUGRA model 

It is easy to see that Aa^ favors light second generation sleptons, while BF(b — > sj) prefers 
heavy third generation squarks. This situation is hard to realize in the mSUGRA model, 
and may be an indication that one must move beyond the assumption of universality, 
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wherein each generation has a common mass at Q = Mqut- In this section, we advocate 
expanding the parameter set of the mSUGRA model to the following values: 



m (l), m (3), m H , m 1/2 , A , tan/?, sign(fi). (3.1) 

In the above, TOo(l) will be a common scalar mass of all first generation scalars at Q = 
Mqut, while mo (3) is the common mass of all third generation scalars at Mqut- tuh is the 
common Higgs mass at Mqut- The above parameter set is well motivated in SO (10) SUSY 
GUT models, where the two MSSM Higgs doublets typically occupy a 10 of 50(10), and 
each generation of scalars, along with a SM gauge singlet N occupies the 16 dimensional 
spinorial representation of 50(10). Thus, the 50(10) gauge symmetry enforces universality 
within each generation, although each generation may be split one from another even when 
50(10) is unbroken. In fact, it is well known that in models where 50(10) is valid beyond 
the GUT scale, at least the third generation of scalars will evolve to independent mass 
values even if all generations begin with a common mass somewhere above Mqut] see for 



instance, Ref. [36]. 



The step of breaking generational universality must be taken with some caution, since 
in general it can lead to violations of constraints from FCNC processes. The most stringent 
of the FCNC limits comes from contributions to the Kl — Ks mass difference, and from 
limits on jj, — ► e'y decay. These bounds, however, apply to splittings between the first and 
second generations. 

Constraints from FCNC processes are usually presented as bounds on off-diagonal 
terms of squark or slepton soft SUSY breaking(SSB) mass matrices or trilinear soft breaking 
terms in the super- CKM basis g |?], H- Here, we are interested in a "proof-of-principle" 



that large splittings between first and third generation scalars can exist and satisfy con- 
straints from FCNC processes. We will assume a scenario of "minimal mixing", where 
SSB mass-squared and trilinear matrices are diagonal at the electroweak scale, and that 
off-diagonal elements are generated only by rotation of the weak scale SSB matrices to the 
super-CKM basis 4 . In addition, in our constraint calculations we neglect the contribution 
of trilinear soft terms, since these contributions are often smaller than those from the mass 
matrices. An explicit case study is shown in the Appendix. 

Within the minimal mixing scenario, the limits from the Kl — Ks = Amx mass 
difference in the super-CKM basis can be translated to limits on mass splittings of diagonal 
soft SUSY breaking mass terms. In the case of equal squark and gluino masses, these 
bounds are roughly [j3§fl 

2 

K(l)-m 9 ~(2)| ~2m c -f , (3.2) 

1Vi w 

which yields mass splittings for weak scale SSB squark masses of order the charm quark 
mass (m c ) if the average squark mass is rriq ~ Myy- If squark masses are much heav- 
ier, the limits become less stringent. Allowed regions of the mg(l) vs. rriq(2) plane are 
shown in Fig. || for several different gluino masses. The limits are obtained from the 
complete gluino-squark box diagram computation, where it is conservatively required that 



4 Arbitrarily large off-diagonal SSB elements would lead to fatal conflicts with FCNC bounds. 
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|Am^ I < Am^ p , and Am^ is the measured Ki — K$ mass difference. The limits are 
presented assuming i.) just a contribution from left squark-left squark mass terms (LL: 
green-shaded region), or more stringently ii.) from the presence of both left-left and right- 
right (RR) squark mass terms (blue-shaded region). The less stringent constraint assumes 
that there is splitting only between soft SUSY breaking masses for left squarks of the 1st 
and 2nd generations, while the corresponding right squarks remain degenerate. In this way, 
non-degeneracy of up to 50% can be allowed for squarks in the TeV range. It can be seen 
from Fig. || that a liberal allowance of non-degeneracy does not work with either scenario 
if squarks ~ 100 GeV, because both areas converge to a line of nearly perfect degeneracy, 
which was the main motivation for the mSUGRA model. In this paper, in order to acco- 
modate the kaon mass difference constraint, we will, without further discussion, maintain 
universality, but just between the first and second generations: mo(l) ~ mo(2). 

Splitting the third generation from the first and second can also potentially lead to 
violations of FCNC processes. One of the main experimentally measured bounds on FCNC 
processes in this case comes from Bjj — B^ mass splitting. As in the kaon case, we neglect 
all but the largest contribution to Amg coming from gluino-squark box diagrams. In this 
case, for equal gluino and squark masses, the soft term mass splitting limit is 



m. 



1)_m9 ~ (3)l Mv? (3 ' 3) 



which is much less restrictive than the kaon case, for both low and high squark masses. We 
adopt the result for gluino box diagrams from Ref. |37]] to calculate the SUSY contribution 
to ArriB d - The result for various values of rrig is presented in the rriq(l) vs. rriq(3) plane in 
Fig. ||, where rriq(i) is the weak scale squark mass for generation i. The magenta regions 
give values of \Am^ SY \ < Am^ p /1000, while green gives \£\m^ JSY \ < Am^^/lOO, and 
yellow gives \£±m^ JSY \ < Am^/lO. From the figure, it is clear that practically the whole 
parameter space of NMH SUGRA model is allowed by the Am^ constraint. 

We have also checked the potential constraint on mo(2) — mo(3) splitting from BF(b — > 
57) measurements. These constraints also yield weak bounds on generational splitting; see 
the Appendix for an explicit calculation. More details can be found, for example, in Ref. [38] 
(see Eqs. (53-54)) and Ref. || (see Eq. (19)), along with discussion of these bounds. 

Independent of which precise limit on Amg we choose, it is easy to see from Fig. |]that 
any constraint on degeneracy of weak scale squark masses applies only if squark masses 
are quite light, well below the TeV scale. As we move to TeV level squark (and gluino) 
masses and beyond, the limits become essentially non-existent. This corresponds to a 
partial decoupling solution to the SUSY flavor problem. 

However, even if squark mass splittings are very large at the scale Q = Mqut, the 
weak scale mass splittings are often much smaller. As an example, we illustrate the running 
of soft SUSY breaking parameters in Fig. ||, where we choose 

m (l) = 100 GeV, m (3) = m H = 1400 GeV, m 1/2 = 550 GeV, 

A = 0, tan^ = 30, fi > 0, m t = 175 GeV. (3.4) 
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• lAm^susv'l < Am^T" for (LL) and [RR) 

• lAm^susvi sAm^pforlLL) only 



Figure 3: Constraints on first and second generation squark masses from Aitik- 

While first generation slepton masses are renormalized to values of ~ 100 — 300 GeV, 
the first generation squark masses evolve to weak scales values in the TeV range, since the 
QCD contributions to RG running are large 5 . Meanwhile, the running of third generation 
squarks is actually suppressed somewhat, owing to the large top quark Yukawa coupling. 
The result is that, even beginning with a large splitting in generations at Q = Mgut, a 
much less severe splitting amongst squarks may be obtained at the weak scale. Of course, 

A simple approximate formula relating weak scale to GUT scale squark and slepton masses is that 
m? ~ mo + (5 — 6)m^ 2 , while m| ~ ttlq + (0.15 — 0.5)m^y 2 . 
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Figure 4: Constraints on first and third generation squark masses from Am^. 



the splitting amongst sleptons remains rather large, and could give large contributions to 
r — * e7 and r — * fij decay. However, the current experimental limits from flavor changing 
radiative r decay are not overly constraining^]. 

Motivated by above considerations, we next generate random points in the parameter 
space given in Eq. 3.1, and again calculate the resultant x 2 values. We have scanned over 



the following region of parameter space for sign(fi) > 0: 

mo(l) : - 3 TeV, m (3) : - 3 TeV, m H : - 3 TeV, 
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m (1,2) = 100GeV, m (3) = 1400GeV, m 1/2 = 550GeV, 
A = 0, tanp = 30, n > 0, m, = 1 75GeV 
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Figure 5: Evolution of soft SUSY breaking masses in the NMH SUGRA model. 

mi/2 : - 1 TeV, Aq : —3 — 3 TeV, tan/3 : 2 - 60. (3.5) 

The colors of the points correspond to those of Fig. [I], so that green gives low x 2 ) red gives 
high x 2 i an d yellow gives intermediate values. A variety of frames showing correlations 
amongst the parameters are shown in Fig. || Only points satisfying the LEP2 constraints 
have been plotted. This is why, for example, mj/2 is bounded from below at ~ 200GeV, 
corresponding to the limit on the chargino mass of > 103.5 GeV. 

From the top row of figures, we see that rather low values of mo(l) ~ — 400 GeV are 
preferred. These parameter values give sufficiently light smuon and mu sneutrino masses 
so as to fulfill the (g — 2) M constraint. Values of my 2 ~ 400 — 800 GeV are preferred, 
although little preference is shown for Aq. One can also see little preference for tan/3 so 
long as mo(l) remains small. For larger values of mo(l), there is a preference for very high 
values of tan (3 ~ 50 — 60, as in the mSUGRA model case. 
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The plots of the second row show that there is some preference for ran ~ mo (3), and 
that there is a preference for mo(l) <C mo (3). The third row shows that while mo(l) ~ 
0—400 GeV, mo (3) ~ 500—3000 GeV is preferred. Further, the last frame of this row shows 
that while low x 2 points can occur at any tan/3 value if there is a large m-o(l) — mo (3) mass 
splitting, that mainly large tan (3 is preferred if the generational mass splitting is small 
(which takes us back towards the mSUGRA case). If may also be pointed out that SUSY 
IMH models are greatly disfavored. Finally, the last row of plots shows the mechanism 
for annihilating neutralinos in the early universe. The first frame shows a region of low 
X 2 where 2m^ ~ hia, i-s. the ^4-funnel at large tan/3, which occurs in the subset of 
mSUGRA like models. The second frame shows that some models get rid of neutralinos by 
stau co-annihilation, while the third frame shows that now many models may also destroy 
neutralinos in the early universe via Z\ — e and Z\ — jl co-annihilation. 

Motivated by the above scan, we next adopt the value tuh = mo (3) (to reduce param- 
eter freedom), and plot the \fx 2 values in the mo (3) vs. m,\/2 plane in Fig. [j], for mo(l) 
values of 50, 100 and 200GeV, with tan/3 = 10, Aq = and /i > 0. The corresponding 
contour plots of BF(b — ► 57), and fm 2 are shown in Fig. ||. 

The most striking feature of the plots is that most of the area displayed is excluded. 
In this case, slepton masses are quite light, in the vicinity of a few hundred GeV. As m^j 
increases, ultimately mg becomes greater than rug, and one violates the cosmological 
constraint on stable charged relics from the Big Bang. Of the surviving parameter space, 
we find significant regions with relatively low x 2 ■ The plot with mo(l) = 50 GeV has a 
rather broad band of low x 2 ■ I n this case, neutralinos in the early universe can annihilate 
by a combination of t-channel slepton exchange (as in the bulk region of mSUGRA), and 
by neutralino-slepton co-annihilation. In addition, smuons and mu sneutrinos are relatively 
light, giving a large, positive contribution to Aa M , while top squarks inhabit the TeV and 
beyond range, effectively suppressing anomalous contributions to BF(b — > S7). In addition 
to the region with low values of mo (3) and low rni/2, a thin strip of the HB/FP region 
survives, where the small value of \x that is generated leads to a Higgsino-like LSP. In this 
case, the HB/FP has moderate but not low x 2 values because the contribution to Aa M is 
low. Finally, we note in this plot that the stau co-annihilation region is somewhat enhanced 
as well on the left most boundary of parameter space. This is because neutralinos can 
annihilate via a combination of t-channel slepton exchange, and also stau co-annihilation. 

If we shift to the mo(l) = 100 GeV plot, then slepton masses correspondingly increase, 
suppressing the relic density in the low mi/2 region of the plot. One must move close to 
the slepton co-annihilation region at the upper boundary of parameter space in order to 
generate a sufficiently low value of Qg^h 2 in accord with the WMAP analysis. Large mo (3) 
values are favored over small ones only because small mo (3) implies larger mj/j values to 
get the relic density right, but then the large value drives up the smuon masses and 
thus suppresses Aa^ too much. The further shift to mo(l) = 200 GeV leads to a plot in 
which very little area contains a low value of x 2 ■ This is due to the fact that slepton masses 
have increased, and neutralino annihilations graphs with t-channel slepton exchange are 
suppressed. Only the border region survives, where p, — Z% and e — Z\ co-annihilation 
occurs at a large rate. 
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Figure 6: Scan over parameter space of NMH SUGRA model. Green points have low x 2 , while 
red points have high \ 2 \ yellow points have itermediate x 2 ■ 



In Fig. [9|, we again show the mo(3) vs. roiw plane, this time for tan/? = 30, and for 
mo(l) values of 50, 100 and 200 GeV. Fig. ^ presents the associated BF(b — » 57), and 
f2/i 2 contour plots for the same tan f3 value. 

By increasing tan/3, we also increase the SUSY contribution to a M for a given set of 
slepton and chargino/neutralino masses. The result is a band of very low \ 2 points in the 
mo(l) = 50 and 100 GeV plots, where essentially a perfect fit to (g — 2) M , BF(b — » 57) 
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Figure 7: Plot of regions of ^px 2 in the mo(3) vs.m 1 / 2 plane of the NMH SUGRA model for 
mo(l, 2) = 50, 100 and 200 GeV with Ao = 0, /i > 0, and tan/? = 10. The green regions have low 
X 2 , while red regions have high \ 2 . Yellow is intermediate. 



and ^z^ 2 can ^ e obtained. The low m^j portion of the mn(3) vs. my 2 plane are largely 
excluded because they give rise to too large a relic density of 0.129 < < 0.5, and 

because they give rise to too large a value of Aa^. In these plots, we also note that the 
HB/FP strip has moved to a somewhat lower x 2 value compared to the tan/? = 10 plots 
due to a larger SUSY contribution to Aa^. Sample spectrum for a point with very low 
X 2 is shown in Table |]. The last frame, with mn(l) = 200 GeV, has become again largely 
excluded, save for a narrow strip where slepton co-annihilation occurs, and in the HB/FP 
region. 
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Figure 8: Contour levels for BF(b — + S7), a p and flh 2 in the mSUGRA model in the mo(3) vs.m\/2 
plane of the NMH SUGRA model for m (l,2) = 50, 100 and 200 GeV with A = 0, fj, > 0, and 
tan/3 = 10. Magenta region excluded by LEP2 searches. 



4. Implications for colliders 

In the previous section, we have seen that the adopton of non- universal scalar masses can 
lead to a reconciliation of the SUSY explanation for a^, BFib — > 57) and 0,^ h 2 . The 
scenario of NMH SUGRA model works if mo(l) — ?rio(2) *C mo (3), and leads to spectra 
typically with squarks and third generation sleptons in the TeV range, while first and 
second generation sleptons have masses in the range of 100-300 GeV. The presence of 
rather light first and second generation sleptons in the sparticle mass spectrum in general 
leads to enhancements in leptonic cross sections from superparticle production at collider 
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Figure 9: Plot of regions of ^fx 2 in the mo(3) vs.m 1 / 2 plane of the NMH SUGRA model for 
mo(l, 2) = 50, 100 and 200 GeV, with Ao = 0, /i > 0, and tan/3 = 30. The green regions have low 
X 2 , while red regions have high x 2 ■ Yellow is intermediate. 



experiments, compared to the case where selectrons and smuons are in the multi-TeV 
range. The enhancement comes from the fact that sleptons may now be produced with 
non- negligible cross sections at colliders, and also in that their presence enhances the 
leptonic branching fractions of charginos and especially neutralinos. In this section, we 
discuss the implications of light selectrons and smuons for the Fermilab Tevatron collider, 
the CERN LHC and a linear e + e~ collider operating with ^/s ~ 0.5 — 1 TeV. 
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Figure 10: Contour levels for BF(b — » 57), a M and Slh 2 in the mSUGRA model in the 
mo (3) vs.mx/2 plane of the NMH SUGRA model for m (l,2) = 50 GeV, and m (l,2) = 100 
GeV with A = 0, fi > 0, and tan/3 = 30. Magenta region excluded by LEP2 searches. 



4.1 Tevatron signals 

At the Fermilab Tevatron collider, the total production cross section for slepton pair pro- 
duction is rather low (~ 10 fb for ~ 100 GeV)p9|], so that prospects for their direct 
detection are not encouraging [4C] even with optimistic projections for the integrated lumi- 
nosity to be gathered. Even so, if sleptons are light enough, then charginos and neutralinos 
may directly decay into them via two body modes: W\ — ► vi&l, and Z2 — > £r£, V£f£ 
and £l£ (we suppress bars on anti-particles). Even if two-body decays are not kinemati- 
cally open, then relatively light sleptons can yield enhancements in chargino and especially 
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neutralino three body decay to leptons. 



If sparticles are accessible to the Fermilab Teva- 
tron, then it is usually expected that pp — > W^W^X 
and W1Z2X will be the dominant production cross sections [4^ 
If W\ — > tvtZ\ and Z2 — ► then clean trilep- 

ton signals may occur at an observable rate [ 42 1. Sig- 



parameter value (GeV) 

~M~2 



nal and background rates have recently been investi- 
gated in Ref. [43], where Tevatron reach plots may 



also be found. Here, we generate all sparticle produc- 
tion processes using Isajet 7.69 for the parameter space 
point m (3) = 1400 GeV, m 1/2 = 225 GeV, A = 0, 
tan (3 = 10 and /i > 0. We plot the isolated trilep- 
ton signal after using cuts SC2 of Ref. Q], where the 
backgrounds are also evaluated. The results, plotted 
versus variation in the mo(l) parameter, are shown in 
Fig. 11, where the signal level needed for a 5a signal 
with 10 fb _1 is also denoted. The error bars show the 
Monte Carlo statistical error. When mo(l) = mo(3) 
(at mo(l) = 1400 GeV), the results correspond to the 
mSUGRA model, and the isolated trilepton signal is 
well below discovery threshold. As mo(l) decreases to 
smaller values, the isolated trilepton rate drops. This 
is due in part to a slight reduction in total SUSY cross 
section (shown in frame &)), but also due to an inter- 
ference in the Z2 leptonic decay rates, due to destruc- 
tive interference between slepton and Z boson mediated 
graphs. As mo(l) drops to even lower values, the light 
sleptons begin to dominate neutralino three body de- 
cays rates, and consequently the trilepton cross section 
rises steeply, to the level of observability. Eventually 
chargino and and neutralino two body decays to slep- 
tons turn on (in this case, first Z2 — > £r£), and trilepton 
rates become very high. For even lower mo(l) values, 
neutralino decays to i>(yi turn on, and briefly suppress 
the trilepton rate, until finally Z2 — ► turns on, and 
the trilepton rate picks up again for the lowest mo(l) 
values. 
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351.1 
184.2 
516.9 
1067.7 
939.8 
910.0 
943.5 
907.1 
1175.1 
1477.5 
1460.0 
1637.1 
319.3 
188.2 
295.1 
1386.1 
1475.4 
1468.5 
348.2 
542.4 
179.4 
347.2 
1379.3 
118.4 
0.115 
3.52 x 10" 4 
35.1 x 10~ 10 



Table 1: Masses and pa- 

rameters in GeV units for 
mo (3), TO1/2, A a , tan/3, sign(fi) = 
1500 GeV, 450 GeV, 0, 30, +1 in 
the NMH SUGRA model. We also 
take run = TOo(3) and TOo(l) = 100 
GeV. The spectrum is obtained 
using ISAJET v7.69. 



The large rates for clean trileptons at the Fermilab 
Tevatron occur for mo(l) values below about 200 GeV, which is just the regime needed 
to give substantial SUSY contributions to (g — 2)^. In the plot shown, the value 
chosen was small enough that sparticle total production cross sections are large enough to 
generate an observable signal. For larger m\/2 values, an observable signal at the Fermilab 
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Figure 11: Rates for isaolated trilepton events at the Fermilab Tevatron pp collider, after cuts 
SC2 from Ref. [El. 



Tevatron doesn't always occur, even if mo(l) is small. 
4.2 CERN LHC 

The reach of the CERN LHC for SUSY particles in the mSUGRA model extends up to 
m i/2 values as high as 1400 GeV (700 GeV) for small (large) values of the parameter mo, 
assuming 100 fb _1 of integrated luminosity [fill. This corresponds to a reach in m g of 3 (1.8) 
TeV, respectively. For the case of the NMH SUGRA model, m^ 2 is bounded from above 
by the requirement that the neutralino, and not a slepton, is the LSP. The upper bound on 
mi/2 is well below the LHC reach in m\ii for the mSUGRA model, and so the CERN LHC 
should easily establish a signal in all of the NMH SUGRA model parameter space with 
low enough mo(l) values as to satisfy the (g — 2)^ deviation. A possible exception occurs 
if SUSY lies in the HB/FP region. But even here, first and second generation sleptons are 



relatively light, and may be accessible to LHC searches in the dilepton channel [46]. 

The collider signals for SUGRA-like models at the CERN LHC are naturally divided 
up according to number of leptons in the final state. Thus, in Ref. ^M, signals for jets 
plus f^T plus 0, 1, 2 same sign (SS) or 2 opposite sign leptons (OS), 3 leptons and 4 or 
more leptons occur. In Fig. we plot signal rates from the NMH SUGRA model for 
mo (3) = 2300 GeV, m 1/2 = 300 GeV, A = 0, tan/3 = 10 and n > 0, versus m (l). The 
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cuts are listed in Ref. |4q| , and have been optimized to give best signal-to-background ratio 
for mo(l) = 100 GeV. For very large mo(l) — mo(3) (mSUGRA model case), a variety 
of nonleptonic and multilepton signals occur, all at observable levels above background 
estimates. As mo(l) decreases, the total sparticle production cross section increases sub- 
stantially, mainly because first and second generation squarks are decreasing in mass, and 
enhancing strongly interacting sparticle pair production rates. As mo(l) decreases further, 
leptonic rates also increase, in part because of increased total production cross sections, but 
also due to enhanced chargino and neutralino (s)leptonic branching fractions. As mo(l) 
drops below about 200 GeV (the value typically needed to explain the (g — 2) M anomaly), 
multilepton rates rise steeply. Thus, we would expect that SUSY as manifested in the 
NMH SUGRA model would be easily discovered, and what's more, the signal events would 
be unusually rich in multilepton events. Such multilepton events can be especially useful 



for reconstructing sparticle masses in gluino and squark cascade decay events [47]. 
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Figure 12: Rates for various SUSY signals at the CERN LHC, using cuts from Ref. |l5). The 5a 
discovery cross section for the $t+ jets channel for 100 fb _1 of integrated luminosity is 53.4 fb. 

4.3 e + e~ linear collider 

A linear collider operating at y/s = 0.5 — 1 TeV may be the next frontier particle physics 
accelerator beyond the CERN LHC. Depending on sparticle masses and the collider energy, 
charginos and neutralinos may or may not be accessible. However, in the NMH SUGRA 
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model, light first and second generation sleptons are needed both to explain the {g — 2) M 
anomaly, but also to enhance neutralino annihilation in the early universe. This means 
slepton masses are typically in the 100-300 GeV range, and likely within reach of a linear 
e + e~ collider [48]. 



5. Conclusions 



In this paper, we have examined the parameter space of gravity mediated SUSY breaking 
models, to see if constraints from (g — 2)^, WMAP Sl^/i 2 and BF(b — ► 57) can simulta- 
neously be satisfied. In performing this task, we used new theoretical and experimental 
determinations of (g — 2)^ which make use of e + e~ — > hadrons data at low energy to 
fix the hadronic vacuum polarization contribution to = (g — 2) M , and which gives a 
~ 3a discrepancy with the SM. While {g — 2)^ favors a spectrum of relatively light sec- 
ond generation sleptons, the BF(b — > 57) measurement favors relatively heavy (TeV scale) 
third generation squarks. Satisfying all constraints is possible in the mSUGRA model in 
only very limited regions: the stau co-annihilation region at low tan/3 ~ 10, and in the 
^-annihilation funnel at very large tan/3. In particular, the HB/FP region gives only a 
tiny contribution to (g — 2) M , and is disfavored by about 3a if the e + e~ data is used for 
the hadronic vacuum polarization determination. 

We advocate relaxing the universality assumption between third generation scalars and 
their first and second generation counterparts. Simultaneously fulfilling all three constraints 
favors low values of mo(l) — ^o(2) ~ 50 — 200 GeV, while mo(3) values of a TeV or beyond 
are preferred. This comprises a normal scalar mass hierarchy, or NMH, at the GUT scale. 
The weak scale mass splittings amongst the squarks can be modest in NMH SUGRA 
scenario, even though mass splittings at Mqut are large. The FCNC constraints from Bj — 
Bd and BF(b — > 57) are very weak and hardly affect the NMH SUGRA model parameter 
space. However, the weak scale splitting amongst sleptons is large, with selectrons and 
smuons in the 100 — 300 GeV range, while staus are in the TeV range. Neutralinos in 
the early universe annihilate via a combination of i-channel slepton exchange and slepton 
co- annihilation. 

A testable consequence of the NMH scenario is that sleptons may be directly observ- 
able at the CERN LHC, provided that the slepton-neutralino mass gap isn't too small. 
In addition, the presence of relatively light sleptons enhances the rates for multilepton 
production in SUSY cascade decay events at the CERN LHC. The light sleptons might 
also lead to observable rates for trilepton events at the Fermilab Tevatron collider. Finally, 
selectrons, smuons and their associated sneutrinos in the 100 — 300 GeV mass range should 
be acccessible to a linear e + e _ collider operating at -y/i = 0.5 — 1 TeV, even if squarks and 
third generation sleptons are beyond the limit for direct searches. 
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Appendix 

In this appendix, we decribe the "minimal mixing" scenario for soft SUSY breaking mass 
terms we adopted for our study, and apply it to constraints from BF(b — * 57) ||. (Our 
notation is similar to that of f38|.) The relevant Lagrangian terms are 



C 3 



~Qi ( m Q)ij 



(5.1) 



where i,j = 1 — 3 are generation indices and a is an SU(2) index. We will ignore the 
possibility of large CP violating phases for this sample analysis. 

Constraints are generally presented in terms of soft SUSY breaking matrix elements 
in the super-CKM basis, where quark mass matrices are diagonalized by 



(5.2) 



and where the CKM matrix is given by K = V£V L . We adopt a basis wherein f u is already 
diagonal, so that V$ = V£ = 1 and K = . 

At Q = Mqut the SSB mass-squared matrix has the form 



1117 



mg(l) 





mg(l) 








m 2 (3) 



(5.3) 



while the trilinear matrix has only one non-zero entry given by (a^)33 = fbAo, where fb 
is the 6-quark Yukawa coupling. 

In the minimal mixing scenario, the generation of off-diagonal terms in the running 
from Mqut to M wea k is neglected. Therefore all off-diagonal elements in the down type 
squark mass matrix sub-blocks {m?~)LL and (m|)^ arise from the rotation to the super- 
CKM basis. These off-diagonal elements are of the order (A^)ll ~ 1CT 2 x (m%)LL while 



(A 



d)LR 



10 



-4 



x (rrijn)LL- For example, for our sample parameter space point given by 



|3.4) these matrices will have the following numerical form 



( m j)iL 



1.56 x 10 b -39.3 2550 
-39.3 1.56 x 10 6 -1.8 x 10 4 
2550 -1.8 x 10 4 2.0 x 10 6 



(5.4) 



and 



{ m d)LR 



-0.1 0.9 -22.2 
0.9 -6.6 160.3 
-22.2 160.3 -3909 



(5.5) 



where all entries are in units of GeV 2 . 

It is also possible to go a step further than the minimal mixing scenario and include 
(two-loop) RG evolution of SSB matrices nig and from Q = Mqut to Q = Mz to 
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generate flavor violating soft terms. We have checked using SoftSUSY[49| that these off- 
diagonal elements would be of the same order of magnitude as off-diagonal elements of SSB 
matrices in the minimal mixing scenario after rotation to super-CKM basis. Therefore, 
for our order-of-magnitude estimation the minimal mixing scenario appears justified. For 
example, for our sample parameter space point ( |3.4| ), the SSB matrices after RG running 
from Mqjjt to Mz will have the following numerical form: 



1.56 x 10 6 -43 -660 
-43 1.56 x 10 6 -7625 
-660 -7626 2.0 x 10 6 



(5.6) 



(units of GeV 2 ), and 



ad 



-2.1 -4.1 -2.3 
-4.1 -20.0 -28.1 
-2.2 -25.9 -665.1 



(5.7) 



(in units of GeV). 

A rotation of the down squark sub-block mass matrices to the super-CKM basis yields 



1.56 x 10 6 -164 3601 
-164 1.56 x 10 6 -2.6 x 10 4 
3601 -2.6 x 10 4 2.0 x 10 6 



(5.8) 



and 



-6.7 0.04 
0.04 -116.4 
-0.66 4.2 



1.1 

-8.4 
-3922.; 



(5.9) 



In this example at least, it is clear that the off-diagonal elements after including RG 
evolution are of the same order of magnitude as those in the minimal mixing scenario. 

One should also notice that (A^ll is proportional to rn^ 3 — m 2 dl , which in our scenario 
is of the order (TeV) 2 , while the {A^lr elements are at most of the order (10 GeV) 2 and 
therefore contributions from A terms can be neglected in our estimates of constraints from 
Arrix and Am#. 

The general bounds from BF(b — > 57) on these matrix elements from 



arc 



,23\ - o n I m g \ 2 2 



{Af) LL < 8.2 



TTi ~ 

500 GeV J q 



(5.10) 



which for (Af ) LL ~ 10~ 2 x m~ implies rriq > 17.5 GeV, which is easily satisfied. Also, 



(A 23 )^ < 0.016 (^§^) 2 ^~ ~ 64 GeV 2 (5.11) 

for rriq ~ 1 TeV. Thus, the assumed mass splittings in the NMH scenario seem safe from 
FCNC constraints when minimal mixing is assumed. 
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